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1. Introduction 

Recent investigations in pEj reduced the study of the nonlinear Schrodinger 
equation (NLS), 

(1.1) iu t + Am — \u\ 2 u = 0, u(0, x) = Uo(x) for x G T 2 , 
to the "Toy Model" dynamical system given by the equation 

(1.2) - id t bj(t) = -\b 3 {t)\%(t) + 26?_ x 67(t) + 2b] + Mt) 
for j = 1, . . . , N, with boundary conditions 

(1.3) b (t) = b N+1 (t) = 0. 

The b/s approximate the energy of families of resonantly interacting 
frequencies to be described below. The main purpose of this paper is 



to study the evolution equation (1.2), both to gain additional insight 



into (1.1) and for its own sake. 



In addition to showing how (1.2) approximates (1.1), a key result of 



[I] is the construction of a solution to (1.2) which transfers mass from 
low index j to high j. In the underlying NLS problem, this implies there 
exist arbitrarily large, but finite, energy cascades. Thus, [I] showed 
that Hamiltonian dispersive equations posed on tori can have "weakly 
turbulent dynamics," the phenomenon by which arbitrarily high index 
Sobolev norms can grow to be arbitrarily large in finite time. 

The question of energy cascades in infinite dimensional dynamical 
systems was considered by Bourgain [2], who asked if there was a so- 



lution to (1.1) with an initial condition uq G H s , s > 1, such that 



^1.4) lim sup ||w(£)||# s = oo. 



t— >OD 



This corresponds to a weakly turbulent dynamic, as there is growth in 



high Sobolev norms, but no finite time singularity. Indeed, since (1.1) 
is defocusing it has a bounded H l norm. One can view this behavior 
as an "infinite-time blowup." 

Although the result in [1] does not answer Bourgain's question, it 
makes significant progress. The result says that given a threshold K ^> 
1 and 5 > there exists uq G H s with ||uo||fp < 5 and T > such that 
||tt(T)||# s > K, where u is the solution to the NLS with u(0) = Uq. 
This establishes 

(1.5) inf < limsup ( sup ||m(;£)||#s J > = oo, 



but not (1.4). This is one of the first rigorous result exhibiting the shift 
of energy from low to high frequencies for a nonlinear Hamiltonian PDE 
viewed as an infinite-dimensional Hamiltonian dynamical system, see 
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also work by Kuksin [9]. The works Carles- Faou [3J, Hani p], and 



Guardia-Kaloshin [6] have also recently treated (1.1). A particular 
achievement of these newer works is their careful construction of error 
estimates on the non-resonant terms. 

The dynamics in [3] were not shown to be generic. Rather, the 
authors constructed a single solution with the desired properties. The 



stability of this solution to the flow (1.2) is unknown. One purpose of 



this note is to explore this question of "genericity" , by investigating 



ensembles of data for (1.2), and finding that, on average, there is a 
transfer of energy from low to high indices. 

In addition to this statistical study, we seek out other interesting 



dynamics in (1.2). Notable behaviors we found include: 

• Compactly supported, time harmonic, structures; 

• Spatially and temporally periodic solutions subject to the adop- 
tion of periodic boundary conditions, 

(1.6) b (t) = b N (t), b N+1 (t) = hit); 

• Nonlinear hyperbolic behavior with both rarefactive waves and 
dispersive shock waves. 

Many of these solutions are obtained by going to the hydrodynamic 
formulation of the problem. Making the Madelung transformation, 

(1.7) b 3 (t) = ^t)exp(ifa(t)) 

with pj > and <pj G R, we obtain evolution equations for pj and <fif 

(1.8a) fa = -pj + 2p j __i cos [2(0j_! - fa)] + 2p j+1 cos [2(fa +1 - fa)] , 
(1.8b) pj = -4pjPj-i sin [2(fa^ - fa)} - kpjpj+i sin [2(fa +1 - fa)} . 

From this perspective, it is clear that phase interactions play a key role 
in the dynamics. 

2. Properties of the Toy Model 



In this section, we briefly review the connection between (1.1) and 



(1.2), and review some important structural properties of (1.2). 



2.1. Relationship to NLS. First, we summarize the argument from 
[3] which relates NLS to the Toy Model. This begins by studying NLS 
in Fourier space, 

u (t,x) = an(t)e in - x+lnl2t . 
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After a choice of gauge eliminating certain trivial interactions, the 
Fourier amplitudes {a n } are seen to evolve according to 

(2.1) idiCl n C-nl^ril ^ ^ Q>n\Q"niQ>nz*~" i 

where 

^4 = \n\\ 2 - \n 2 \ 2 + |n 3 | 2 - |n| 2 , 
r(") = {(^1,^2,^3) e (Z 2 ) 3 |ni - n 2 + n 3 = n, ni ^ n, n 3 ^ n] . 

For any n, the most significant contributions in the summation will be 
the elements of T(n) belonging to the resonant set, 

r re sO) = {(ni,7i 2 ,n 3 ) e r(n) | |ni| 2 - |n 2 | 2 + \n 3 \ 2 - \n\ 2 = 0} . 

Restricting ( |2.1 ) to the resonant modes, we have 

(2.2) -id t r n = -r n \r n \ 2 + 

ni,n 2 ,n 3 er iea (n) 



ly ly ty 

' m 1 712 ' 713" 



A union of disjoint sets, Aj, of resonantly interacting frequencies is 
constructed, 

A = Ai U A 2 U ■ ■ ■ U A N , 

where the mass from modes in generation Aj, r ni and r„ 3 , can mix 
to transfer mass to modes r n and r n2 in generation A J+ i, where again 
ni — n 2 + n 3 = n. Subject to certain additional conditions, we will have 
that for all t and j, 

r n {t) = r n >(t), V n,n' G Aj. 

Once these sets have been constructed, a nontrivial step, the relation- 
ship between the toy model and (2.2) is 



(2.3) bj{t)=r n (t), ViigA 



3 ■ 

Hence, |frj(i)| 2 is a measure of the spectral energy density of generation 
A,. 

To show that NLS has an energy cascade, the authors used ideas 
inspired from studies of Arnold diffusion (see [1]) and explicit ODE 
manipulations to show that the Toy Model admits an instability mech- 
anism transferring mass from a low index node to a high index node. 
By the construction of the initial data set A, such a mass transfer in 
the Toy Model yielded growth of high Sobolev norms of the solution 
to the resonant system (2.2), which in turn implied an energy cascade 
for NLS via a stationary phase argument. 
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2.2. Structural Properties. We recall some of the results from Sec- 
tion 3 of [I]. The toy model is a Hamiltonian dynamical system with 
Hamiltonian given by 



(2.4) 



H[b = (h, b 2 ,..., b N )} = ( jh\ A - Mbp. 

7 = 1 ^ 



2 - 



and symplectic structure, 



(2.5) 



.dbj 
dt 



db< 



1,...N. 



This structure applies to both the original Dirichlet boundary condi- 
tions, (1.3), and the periodic boundary conditions, (1.6), studied below. 

The Toy Model, (1.2), admits many of the symmetries of (1.1), in- 
cluding phase invariance, scaling, time translation and time reversal. 
However, many of these symmetries are redundant, and the only known 
invariant, other than (2.4), is the mass quantity, 

N 

(2.6) M[b] = £>| 2 . 

These invariants are useful in assessing the performance of our nu- 
merical schemes. A robust algorithm should preserve them within a 
controllable error. 

Since (1.2) is a finite-dimensional Hamiltonian system, the behavior 
can studied statistically. By Liouville's theorem, the Lebesgue measure 

N N 

j J dbj = Y\ d Re bjd Im bj 



on R 2N is invariant under the dynamics of (2). Moreover, in view of 



the mass conservation, the white noise 

TV 

d^ N = Z- N l e-^^\ 2 \[db ] 



N 



( 2tt) - N Y[ e~ 3 ( Re b > 1 ) 2+ ( Im h > ? d Re bjd Im bj 



is an invariant probability measure for (|1.2|). In particular, the Poincare 

ensures that almost 



recurrence theorem (see, for example, p. 106 of 
every point b in the phase space is Poisson stable. That is to say, there 
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exists {t n }^ =1 tending to oo (and another sequence tending to — oo) 
such that 

lim b(t n ) = b 



where h(t) is the solution to (1.2) with b(0) 



Here, "almost every" is with respect to both the white noise fi and 



the Lebesgue measure Yl^i^bj, since they are absolutely continuous 
with respect to one another. Of course, this is only an "almost ev- 
ery" statement, but it does says that the solution of the toy model 
constructed in [3] is destined to return to the original configuration. 

3. Random Phase Interactions and Ensemble Dynamics 

In this section, we present the results of running an ensemble of 
initial conditions. The statistics of the results indicate that there is 
generic movement of mass from low to high nodes. 

For our first ensemble, we took as the initial conditions, 

[VI - e 2 exp{i0j} , j = j*. 

e G (0, 1), 1 < < N and 9j are identical independently distributed 
random variables, 8j ~ £7(0, 27r), the uniform distribution on [0,27r). 



Thus, (3.1) has mass one, with the majority of the mass concentrated 
at j+, and random phases on each node. 

To study the spreading of energy in this system, we introduce the 
Sobolev type norms, h s , defined as 

N 

(3.2) \\Hl = Y,^\ 2 - 

We note that this norm will measure shift of mass to higher j indices, 
but is difficult to connect directly to the corresponding H s norm of a so- 



lution to (1.1) on the torus since that requires specifying the placement 
function from [I]. 
Another candidate is 

N 

(3.3) (^2^|6,| 2 )^ 

3=1 

which can be more closely connected to the construction in [I]. How- 
ever, since we only simulate a finite number of generations, if one grows, 



the other must too. Thus, we employ (3.2). 



We now proceed with our simulations for (3.1). Taking e = .1 



iV = 100 and = 10, we simulated this initial condition until t = 10000 
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2000 4000 6000 8000 10000 

t 

(a) a = 1 



2000 4000 6000 8000 10000 



(b) a = 2 




(c) a = 3 



(d) s = 4 



Figure 1 . Here we show ensemble averages for the evo- 
lution of normalized Sobolev norms for the initial con- 
ditions (O) with e = .1, JV = 100 and = 10. The 



ensemble size was 10000 and the dashed lines correspond 
to 95% confidence intervals. 



with 10000 realizations of the random phases. Generic slow growth in 
Sobolev norms appears in Figure [T] These were computed using the 
explicit Runge-Kutta Prince-Dormand (8, 9) method, with a relative 
error tolerance of 10" 12 and an absolute tolerance of 10~ 14 , see [5]. 
Over the entire ensemble, the maximum absolute and relative error in 
the invariants remained below 10~ 9 . In Figure [2j we show the evolution 
of a particular realization to show how this growth in norms occurs. 
As the figure shows, there is a spreading of the mass away from the 
initial site of high mass. Additionally, there is local exchange between 
sites. 
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Figure 2. The evolution of a single realization from the 



ensemble with initial condition (3.1). 



Of course, as N — > oo, the nodes in (3.1 ) decouple. As an alternative, 
we consider 



(3.4) 



fexp-fiflj-}, 

, exp{i9j}, j = j+. 



The decay in j ensures that as — > oo, the solution has finite mass. 



The results, plotted in Figure 3 are similar to the (3.1) case. There is 



somewhat more growth in h 1 , but less growth in the other h s norms. 



4. Particular Solutions 



In this section, we consider several particular solutions to (1.2), in- 
cluding localized solutions, periodic solutions and a "hyperbolic" solu- 
tion. These were motivated by the hydrodynamic formulation of the 



problem, (1.8). 



4.1. Localized, Uniform Phase Solutions. We first seek solutions 
which stay in phase for all time, 



(4.1) 



<t>j{t) = <Pj+i{t). 



Such solutions are said to be phase locked. Assuming this holds, (1.8) 
becomes 



(4.2a) 
(4.2b) 



Pj = 0- 
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t 

(a) a = 1 
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8000 10000 



(c) a = 3 



(b) a = 2 




(d) s = 4 



Figure 3. Plots of normalized Sobolev norms for the 
initial condition (3.4) with e = .1, N = 100 and = 



10. The ensemble size was 10000 and the dashed lines 
correspond to 95% confidence intervals. 

We now need a solution to 

(4.3) - P j + 2pj_i + 2p J+1 = loeR, for j = 1, ... N, 

where u is independent of t and p = Pn+i = 0. This corresponds to 
the linear system 

(4.4) 



/- 1 

2 




2 







\ o o 



LOl 



•1/ 



Though the matrix is tri-diagonal, it is not diagonally dominant, so its 
solvability is not immediately clear. However, 



Theorem 4.1. The matrix in (4.4) has no kernel for any N. 
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Proof. Letting An denote this matrix, proving it has no kernel is equiv- 
alent to showing det An 7^ for any N. Indeed, 



(4.5) 



det A 



N 



•1 • det A N -i - 2 • 



2 2 
-1 
2 





2 





— det An-i — 4 det An-2, 



■■ 

... 

... 

2 -1 



giving us a recursion relation for the determinant. By inspection, 
(4.6) det Ai = -1, det A 2 = -3. 



By induction, recursion relation (4.5) demands that for all N, the de- 
terminant of An must be odd. Hence, it is never zero. □ 

A consequence of this is the following Corollary. 

Corollary 4.1. Any nontrivial phase locked solution has u 7^ 0. More- 
over, given any real valued function u(t) 7^ for all t 6 R ; there exists 



a nontrivial phase locked solution {(<j)j , Pj)}j=i to (1.8). 



For a given N, the linear system can always be solved, and a phase 



matched solution of (1.8) exists. However, this will not always yield 
a solution of (1.2). As Figure 111 shows, at N = 5, the solution is not 



strictly positive, and the Madelung transformation cannot be inverted. 
Despite the obstacle at N = 5, we can again obtain a strictly positive 
solution at N = 8 and higher, as Figure [5] shows. It is possible to ask 
whether or not, there exist positive solutions for specific, but arbitrarily 
large, values of 

Since the equations for the toy model, the hydrodynamic formulation 
and (4.4) are autonomous in j, one can concatenate these localized 
solutions together to form new solutions. For example, one could place 
the N = 2 solution at lattice sites 55 and 56, and the rest to zero, 
within (1.2) with iV = 100. Thus, one can construct explicit solutions 



with isolated regions of support on arbitrarily large systems. 



Indeed, using recursive linear algebra techniques, this question has been an- 
swered affirmatively thanks to an observation of Stefan Steinerberger about the 
recursive nature of the sequence of N for which positive solutions exist and 
a. clever argument from Sergei Ivanov through the Math Overflow Project at 
http : //mathoverf low . net /quest ions/ 106816 
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Figure 4. Solutions of (4.4) with to = 1. The solution 



at ./V = 5 is not strictly positive. 
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(b) N = 142 



Figure 5. Solutions of ( |4.4[ ) with u = 1. Both solutions 
are strictly positive 
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4.2. Time Harmonic Periodic Solutions. Here, we consider the 



(4.8a) 



problem with periodic boundary conditions, (1.6). Assume for all j 
the initial condition satisfies 

(4.7a) p i+1 (0) =/>,-_!(()), 

(4.7b) i+1 (O) = i _ 1 (O). 

Furthermore, assume one or both of the variables, p or <fi, are not 
uniformly constant. This will result in time harmonic solutions. This 
can be observed by computing 

= ~ (fij+i ~ Pj-i) 
+ 2pj [cos(2(0 i+1 - (pj)) - cos(2(0 i - (pj-i))} 

+ 2p j+2 cos(2(0 j+2 - <Pj+l)) 

- 2pj- 2 cos(2(0 j _i - 4>j- 2 )), 

j t (Pj+i ~ Pj-i) = ^PjPj+i sin(2(0 i+ i - (pj)) 

(4.8b) ~ 4 Pi+2Pj+i sin(2(0 i+2 - <j> j+ i)) 

- Apj-xpj-2 sin(2(0 i _i - <pj-2)) 
+ 4p i p j „ 1 sin(2(^ - <f>j-i)). 

If, initially, (pj+\ = <Pj-i and pj + \ = pj-i for all j, then they will remain 
so. We also have: 

• pj + Pj-i is constant in j and t; 

• cos(2(0j — <f>j-i)) is constant in j, but may vary in t; 

• sin(2(0j — 4>j-i)) is constant in j, but may vary in t; 

• <pj + 4>j-\ is constant in j, but may vary in t. 

These symmetries reduce the problem to four variables, pj, Pj-i, (pj, 
and 4>j-i- 

(4.9a) ^ = -pj + Apj_ x cos(2(0 i - 0j_i)), 

(4.9b) = - Pj _ x + A P j cos{2((pj - (pj-i)), 

(4.9c) pj = SpjPj-i sin(2(0 i - (pj-x)), 

(4.9d) pj-i = SpjPj-! sin(2(0j - (pj-i)). 
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Defining 






y±. rua ) 




<P = <Pj + 9j 


(4.10b) 




Acp = (pj — (pj 


(4.10c) 




p = Pj + pj. 


(4.10d) 




Ap = pj — pj^ 

r r J r J 


we have: 






(4.11a) 


It* 


= —p + 4pcos( 


(4.11b) 




= -Ap - 4Ap 


(4.11c) 


d 

dt P 


= 0, 


(4.11d) 


d a 
dt 


= 16pjPj-i sin( 



-1- 



Since p 2 — Ap 2 = Apjpj_i and p is invariant, we have a closed system 
of 2 equations for A</> and Ap. 

(4.12a) A(p = -Ap [1 + 4 cos(2A</>)] , 

(4.12b) -^Ap = 4(p 2 - Ap 2 ) sin(2A0). 

(JjL 

The system is Hamiltonian with 

(4.13) H= ^(l + 4cos(2A0)) (p 2 - Ap 2 ) 
and symplectic structure 

(4.14) — A(p = — — , —Ap 



dt T <9Ap' dt r dA(p 

Thus, we anticipate time harmonic motion. An example appears in 
Figure |6j where the initial condition is 

(4.15a) pi = p 2 = 1, 

(4.15b) 0! = |, 2 = 0. 

5. Discrete Rarefaction Waves and Weak Turbulence 



In this section, we explore the dynamics when the initial configura- 
tion is given by the out of phase initial condition <pj = (pj + i — |. If this 
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Figure 6. 



phase relation were to somehow persist, the resulting equations hints 
at the discrete Burger's formulation of the hydrodynamic equations 



(5.1) 



0, = 0, 

Pi = -*PjPj-i + ^PjPj+i = -8pj ( 



This has discrete rarefaction and shock wave dynamics. We call (5.1) 
the discrete Burger's equation since, were we to discretize 



Pt 



-8pVp, 



in space and take the gridpoint spacing paramter equal to one, we would 
recover the above equation. We note here that the rarefaction waves 
we observe have similar dynamics to those found in Fermi-Pasta-Ulam 
chains, see [8]. 



However in (1.8), there are additional terms which prevent this phase 



relation from persisting. We show here how the solutions evolve, with 
the initial condition 



(5.2) 



expz{(j - l)7r/4} 



As a first example, we solve (1.2) with the initial condition (5.2) over 
iV = 100 lattice sites. The results appear in Figures [7] and |8j As can be 
seen, the h s norms eventually cease to be monotonic. To see persistent 
growth in the norms, we can look at a system with N = 5000 sites and 



for longer a time; see Figures [9] and 10 

The simulation on N = 5000 reveals that the rarefaction portion of 
the solution has more structure than is apparent in the case of iV = 



100. As shown in Figure 11, the rarefaction wave evolves with several 
different slopes. 
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Figure 7. A combination of discrete rarefaction and 
dispersive shocking. As an initial condition, pj = 1 at all 
values of j. 




FIGURE 8. Growth in the h s norms for the dynamics of 
Figure [7| 
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Figure 9. A combination of discrete rarefaction and 
dispersive shocking. As an initial condition, pj = 1 at all 
values of j. 




FIGURE 10. Growth in the h s norms for the dynamics 
of Figure [9j 



Unfortunately, as N — > oo, (5.2) will not correspond to a finite mass 



solution. Thus, we studied the weighted initial condition, 



(5.3) 



expz{(j - l)7r/4} /j. 
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This, too, results in energy transfer, though it is not monotonic. Several 



frames from this simulation appear in Figure 12 , and the growth of the 
norms can be seen in Figure [13] The norm growth is quite pronounced; 
this may be due to an inability of mass to propagate backwards, against 
the weight 1/j. 

In similar calculations with 

(5.4) b 3 =expi{(j-l)n/4}/(w(3)y. 

for < a < 1, where w(j) — > oo as j — > oo, we observe that some form 
of the rarefaction front propagates forward even with a tail in the higher 
nodes and that the structure of the rarefaction wave persists longer as 
a — > 0. Since the rarefaction wave and the backward dispersive shock 
travel at finite speeds in the simulation, we observe motion to large 



j on much longer time scales for large N. Here, the weights in (5.4) 



allow us to study rarefaction waves in a setting with h 1 norms of order 
1 instead of order N. In addition, we observe that the rarefaction wave 
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Figure 12. Evolution of the initial condition (|5.3|) 



solution is robust even for initial data of the form (5.4), which has less 



back scattering thanks to the smaller jump at the right endpoint. As 
the rarefaction front enters the decaying tail, it does however begin to 
lose some mass at regular intervals, but continues to propagate weakly. 
All rarefaction wave solutions presented in this section result in norm 



growth when mapped back to solutions of ( 1.1 ). This is due to Propo- 
sition 2.1 of [3], which states that a mass shift to the higher nodes 
in the Toy Model results in growth of higher Sobolev norms of the 
corresponding solution to (1.1). This is fundamental to showing the 



importance of tracking the rarefaction wave front moving toward large 
j in (1.2). However, a more detailed result relating to the frequency 



scales at each generation Aj and a better categorization of families of 
resonant frequencies would be required to address this issue in its en- 
tirety and observe a rarefaction front in the resonant frequencies of the 
torus. It is unclear at this point how to directly translate solutions 



with frequency cascades in (1.2) to computationally observable solu- 
tions with frequency cascades leading to H s norm growth for s > 1 in 



(1.1 ). Here, however, we have demonstrated the robustness of solutions 



that move mass in (1.2) from low to high j. 
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FIGURE 13. Growth in the h s norms for the dynamics 



of Figure 12 



6. Continuum Limit & Compacton Type Solutions 

As a final observation, let us introduce the parameter < /i < 1, 
such that 

(6.1) B(xj,t) = bj(t), Xj = hj. 
Taylor expanding, 

(6.2) - id t B = [3B 2 + Ah 2 ({d x Bf + Bd 2 x B) + 0(/i 4 )] B. 
Neglecting 0(/i 4 ) terms, 

(6.3) - id t B = 3\B\ 2 B + Ah 2 Bd x (Bd x B) . 

This retains the toy model scaling that if B(x,t) is a solution, then 
so is XB(x, X 2 t). It is also invariant to multiplication by an arbitrary 
phase. 



The equation (6.3) is, formally, degenerately dispersive, and it can be 



viewed as an NLS analog of the compacton equations, [131 QUI [HI fT2"] . 
One of the interesting features of the compacton equations is that they 
admit compactly supported nonlinear bound states, which we now seek 



for (6.3). We begin with the ansatz B = e lt Q(x), Q > 0. Consequently, 
Q solves 

(6.4) Q = 3Q 3 + Ah 2 Q(QQ'y 
which can be expressed as 

(6.5) Q = 3Q 3 + 2h 2 Q(Q 2 )" 
Letting U = Q 2 , 

(6.6) 



2h 2 U" + 3U - 1 = 



This can be integrated up once to 

3 



(6.7) 



h 2 U' 2 + -U 2 -u 



c 
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FIGURE 14. The compacton solution given by (6.11) 



h can easily be scaled out by changing the dependent variable, thus we 
set h — 1. We always have a potential well in this equation, so there 
should be homoclinic orbits. 

If C = 0, then < U < 2/3. The explicit solution is 



(6i 



Putting h in, 



Q(x) 



(6.9) 



Qh{x) 



2 . 1 3 
sin —\—x 

3 I 2 V 2 



Phase shifting the solution by 7r/2, we can alternatively have 
(6.10) Qh(^) = y|cosU- 1 iy|x 

We can turn this into a compact structure if we now define 
y^fcos (h-^y/la 



(6.11) Q c h (x) 



\x) \x\ < hm\j\ 



\x\ > hir 



This structure, with h = 1, appears in Figure 14 
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Note that this will satisfy the equation in the strong sense. If we go 
back to the sine formulation and look at x = with h = 1, then 



Q c {x) 



| sin fjy 0<x<nJl 



otherwise. 



Then, near x — 0, 

(6.12a) Q c (x) ~ 

(6.12b) (Q C W) 2 ~ x 2 H(x) 

(6.12c) ^ (gC(x))2 ~ 

(6.12d) Q c (x)-^(Q c (x)) 2 ~ 

where H(x) is the Heaviside function. Hence, the most degenerate term 



in (6.5) is continuous. 
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